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■ Encouraged by experiments on *He in aerogels, we confine planar spins in the pores of simulated 



aerogels (DLCA) in order to study the effect of quenched disorder on the critical behavior of the 
three dimensional XY model. Monte Carlo simulations and finite size scaling are used to determine 
critical couplings Kc and exponents. In agreement with experiments, clear evidence of change in 
the thermal critical exponents f and a is found at non-zero volume fractions of impurities. These 
changes are explained in terms of hidden long-range correlations within disorder distributions. 

PACS numbers: 05.70.Jk, 64.60.Cn, 75.40.Mg 

^ , The superfluid (A) transition of ^He in porous media has been subject of exciting work due to the change of critical 
exponents observed in many experiments Close to the transition temperature (T\), the relative superfluid 

density scales as ps/p^ \T — T\\'' where C ~ 0.671 is the critical exponent for the pure system Ij]. The correlation 
' length behaves as ^~ |T — T^l^"^ where the exponent v is given by the hyperscaling relation |^ I3 

^: c = (rf-2)^, (1) 

d . with d being the dimension of the system. Thus, C — ^ for bulk "^le at the A-transition. Experiments on ^He embeded 
^ in porous media such as silica aerogels and xerogels report an increase of C up to 0.81 and 0.89, respectively 
' Similar studies with other kind of materials, e.g., vycor and porous gold 0jl3|i do not reveal any change of criticality. 
It has been argued that discrepant results between the two families of materials obey different regimes in spatial 



correlations. 



' Harris criterion establishes that short-range correlated (SRC) impurities do not affect criticality in systems with 
a negative specific heat exponent a. The extended criterion due to Weinrib and Halperin |^] (WH) states that disorder 
' can be relevant when it exhibits long-range correlations (LRC), even in the a < case. Quenched defects obeying 
^ I power-law correlations, g{r) ~ at long distances r, are relevant under the condition av — 2 < for a < d (LRC 
condition) where is the exponent for the original pure system. Prudnikov et al. lOj obtained more precise results 
for thermal [v) and magnetic (r/) exponents of general A^-vector models in (i = 3. One of their main results is that 
WH's heuristic prediction, i.e., that z^iijc' = 2/a would be the exponent for the impure system, does not hold in their 
. more accurate field-theoretical description. In practice, the A-transition of '*IIe exhibits a ~— 0.013 4j, so it has been 
' considered as an ideal system to test these criteria 0- 

■ Moon and Girvin 11] (MG) have notably confirmed WH predictions, performing Monte Carlo simulations of the 
three-dimensional XY model (3DXY), which is in the same universality class as the superfluid bulk "^He. These 
authors studied the 3DXY model in the presence of correlated disorder, the fractal incipient percolating cluster, and 
found C — 0.722 in contrast to the pure 3DXY value. They also studied the 3DXY model with random bond dilution 
(which is SRC), and observed no changes in criticality, consistent with Harris criterion. 
(~| Porous materials used in "^He experiments where changes of criticality occur are obtained through sol-gel processes. 

O In principle, they exhibit fractal features only at finite scales (lOnm < r < A ~ lOOnm) and low concentrations 
^ (porosities ip > 90%), appearing homogeneous at large scales or large concentrations According to theory, no 

J> . change of universality should be expected for the superfluid transition in the latter cases: after the WH criterion, as 
T — > Tc (^ ^ 00), disorder should not be relevant once ^ surpasses A. It would seem homogeneous (thus SRC) for 
^ the typical resolution of the system at criticality. Surprisingly, this expected crossover from the LRC regime to the 
d ' SRC regime, where Harris criterion should apply, has never been observed Experimental works reveal 

that these changes of criticality also occur in aerogels with relatively low porosities {(p < 90%) which do not seem to 
be fractal 0|. These results appear to contradict WH criterion, unless there exists LRC within the whole structure 
which is not observable with standard techniques. Furthermore, systematic violations of Josephson's hyperscaling 
relation 0,0 

di/ = 2 - a (2) 

have been evidenced in experiments involving estimations of a and v within the same aerogel sample 0, Q . These 
intriguing results have caused controversy. 
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In this letter, we report Monte Carlo simulations of the 3DXY model confined in the pores of computer modeled 
aerogels. These media are simulated through the on-lattice diffusion limited cluster-cluster aggregation (DLCA) 
algorithm , which has proved to reproduce well the structural features of aerogels . We find a change of critical 
exponents v and a, even at a relatively large concentration, consistent with experiments on '*He confined in aerogels 
P, 0, 13 ■ We argue that hidden LRC within DLCA aerogels would explain these changes in criticality, so results do 
not really contradict the WH criterion. Scaling consistency in experiments and simulations with respect to theoretical 
predictions is established. 

Simulated aerogels are built on simple cubic lattices of different sizes (L = 16 — 64) using a standard DLCA model 
with periodical boundary conditions (PBC) [l^ . A fraction c of the total number of sites in the lattice are 
initially occupied with cL^ randomly disposed particles. These are then allowed to, respectively, undergo Brownian 
motions, until they irreversibly stick together when they enter in contact. Aggregates are also allowed to perform 
Brownian motion, and the procedure is repeated until all monomers and aggregates form a single object. This impurity 
distribution remains quenched in the box in all that follows. 

Empty sites of the lattice (aerogel pores) are filled up with planar spins 4>i = {cos9i,sin6i). No interactions between 
spins and aerogel sites are considered. Ferromagnetic interactions between nearest-neighbor pairs (i, j) are described 
by the XY Hamiltonian with PBC: 

^KH = ■ 0j , (3) 

where K ~ 1/T is the coupling (inverse temperature, with ks = 1); Kij = when the pair includes at least 

one aerogel site; and Kij — K elsewhere. The magnetization squared (order parameter) is given by = 11(0) IP = 
II J2i'j=i ||^/-^^j where N = {l — c)L^ is the effective number of XF-occupied sites. Helicity modulus [(T^)]l^c 
has been determined using the Kubo formula 0, ^3| and averaged over the three directions /t = x, y, z to improve 
statistics, assuming that the model and disorder distributions are isotropic; square brackets express the average over 
disorder realizations for fixed L and c, as will be stated below |l5|. Suitable statistics over disorder were considered: 
from about 256 realizations for L = 16 to 16 for L — 64. A cluster update has been used to minimize the effects 
of critical slowing down |l6l |. defining one sweep as eight consecutive cluster flips. Thermalization is reached after 
2 X 10^ sweeps while 2 x 10^ production sweeps are taken for further statistical analyses. Reweighting techniques 
and finite size scahng (FSS) [H have been performed to obtain critical exponents. Disorder averaging has been 
performed following Ref.|l5l|. Critical couplings Kc have been determined by the crossing of Binder's magnetization 
fourth cumulant, [C/b]l c = 1 — c/3[(A/^)]| ^ , which is size independent at Tc for continuous phase transitions 

Illlig. 

Normalized critical temperatures are shown with circles in Fig.l as a function of impurity concentration c. The 
dashed curve corresponds to the mean field approximation (MFAE), Tc{c)/Tc{0) — (z(c))/6, calculated through the 
average coordination number {z(c)) of the XY sites. This expression approaches our results well at weak values of c. 
Evident deviation at c = 0.2 reflects that higher order terms, beyond the mean field approach, should be taken in this 
case. The random dilution expression (MFRD), Tc{c)/Tc{0) = l—c, is illustrated with the dotted line for comparison. 
The higher coordination number within aerogel pores provide a rough geometrical explanation for the curve concavity 

Cross averages of M" with energy (n=l, 2, 4) are taken to calculate logarithmic derivatives with respect to the cou- 
pling |91n(M")/5if| [13. The correlation length exponents f have been determined using FSS, \dln(^M^)/dK\max^ 
L^/'^ . Power law fits for n = 2 are depicted in Fig. 2 (a) for the pure system and for c = 0.05, 0.1, 0.2. Similar fits have 
been made for n= 1, 4 in all cases obtaining self-consistent results. In agreement with experiments 0,|3,l3i exponents 
in confinement 0.76(3)] differ satisfactorily from the value obtained here for the pure system [iv = 0.67(1)] which 
is consistent with previous works Fits have been performed including data for L>Lmin, varying Lmm from 16 
to 32 observing no change in exponents (within error bars). This reveals that there is no finite size crossover in the 
length scale studied, and suggests that exponents of the 3DXY model confined into DLCA aggregates effectively flow 
to a new universality class. Exponent v decreases slowly in the studied region c = 0.05 — 0.2, though detailed studies 
in denser aggregates should be made to confirm this tendency. 

On the other hand, inflexion points from independent Tl{K) extrapolated curves have been taken to determine 
averages at pseudocritical couplings. Plots for helicity modulus in Fig. 2 (b) ([{T p_)]l^c'^ L~'"^'^) show that the finite 
size scaling exponent v/vc^ 1 is not affected by confinement. Since helicity modulus is proportional to the superfiuid 
density {( — v) 6], the hyperscaling relation |^ is satisfied for d = 3 for all volume fractions tested. Results suggest 
that the effective dimension of aerogel pores occupied by XY spins is Euclidean in fact. 

Specific heat Cl has been computed from energy fluctuations. According to the scaling ansatz ISj, it behaves as 
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(jmax _(jmax ^ j^a/u ^ rpj_^g method developed by Schultka and Manousakis has been implemented to determine 
FSS exponents a/iy directly. While we obtain a/i^ = —0.019(2) [a = —0.013(2)] for the pure system, an approximately 
constant value ajv^ —0.67 yields a~ —0.5 for the confined system, in agreement with experiments on denser aerogels. 
Writing down the relation ||2Jl as i^* = 2/((i+a/z^) we can compare both estimates of the correlation length exponent. 
As shown in Fig. 3, results with v and v* differ in all c ^ cases. This qualitatively agrees with experimental data 
Q IE Q depicted in the figure for comparison. The violation of hyperscaling relation ^ observed in experiments 
is reproduced here in terms of v and v* . However, these results may require a more detailed analysis of finite-size 
effects on the ajv exponent as in Ref . [j^] . 

Despite the agreement between our results and experimental data, an apparent contradiction with the WH criterion 
remains to be explained. As stated above, only LRC within disorder distribution could account for changes of 
exponents observed in simulations and experiments After detailed analyses on aerogels and DLCA aggregates, 

only a finite region of power-law scaling is found for these objects Nevertheless, previous studies show the 

existence of such LRC structures within two-dimensional DLCA aggregates jl^ . We carry out here a detailed study, 
in order to evidence LRC structures within DLCA aggregates in (i = 3, showing that non-fractality of aerogels at long 
scales is only apparent. 

First, we briefly define the different objects studied within the whole structure. The gelling cluster (GC) is the 
first aggregate to reach opposite sides of the box in any direction, tg is the time step at which this occurs; islands 
are then defined as the smaller aggregates not yet attached to the GC at tg. The whole DLCA is finally completed 
just a few steps after tg when islands reach the GC. Additionally, corresponding backbones can be defined for the 
GC (BBGC) and for the whole DLCA (BBDLCA) as performed elsewhere for the two-dimensional case jl^. Using a 
box counting algorithm jl^, we obtain the average density n(r) =m(r)/r'' for r<L/2. It is related to the two-point 
correlation function g{r) ~ r"" as n(r) ~ r~'--'^~'^!\ dj being the fractal dimension of the object and a = 2{d — df). 
Plots of n{r) versus r are sketched in Fig. 4 for all structures described above at volume fractions c = 0.05 and 0.2. 

As observed in Fig. 4, LRC of gelling clusters (squares) appear hidden when one considers the entire distribution 
of impurities at tg, i.e. GC+islands (diamonds) whose plots coincide with those of whole DLCA (triangles). The 
GC alone and respective BBGC (circles) scale at distances comparable to the lattice size {L/2 at PBC). The average 
inverse slope of plots from Fig. 2 (a) has been depicted for a qualitative comparison (dashed line). Arrows indicate 
the crossover to homogeneous regimes, evident for the whole DLCA and the GC-|-islands distributions, while fractal 
regimes of GC and BBGC span all scales. In conclusion from these results, the short-ranged fractal regime of DLCA 
would prevent the structure as a whole to be relevant for criticality of the 3DXY model. 

Authors in Refs.P,|3| noticed the same paradox for the A transition of ""He in aerogels: standard measures evidence 
a finite cutoff for the fractal regime fl2| beyond which aerogels appear homogeneous. If disorder were strictly homo- 
geneous at large scales, then the tail of its correlation function would decay exponentially whence Harris criterion 
applies: a crossover to the SRC regime would be expected as T — > Tc. However, such a crossover is never observed: 
exponents rest still different from those of pure ''He as Tc is closely approached in experiments, and this is due to LRC 
according to the WH criterion. Note that islands have a homogeneous distribution at tg, so when they finally reach 
the GC at random sites, the whole DLCA appears homogeneous in spite of the underlying LRC fractal within. If we 
take the 3DXY model in the presence of a random distribution of impurities as the original critical system, it has the 
same exponents as the pure system. Then, it should certainly be affected by a LRC distribution of vacancies (GC) 
which actually defines the observed critical behavior. WH considered the correlation function as a sum of power-law 
terms, g{r) ^^^Cir~"'' 0, to show that it is the dominant term of g{r) (smallest Oi) which defines criticality. This 
decomposition of the correlation function, together with the analysis of structures sketched above, confirms that the 
change of criticality observed in the 3DXY model is not in contradiction with the WH criterion. 

To conclude, we have found clear changes of criticality in the 3DXY model when confined in DLCA aerogels at 
different porosities. Evidence of long-range correlated structures within DLCA aggregates in three dimensions has 
been presented, applicable to real silica aerogels as well. This correlation functions analysis should clarify the reason 
these changes of criticality occur in the 3DXY model in DLCA, and why the crossover to bulklike behavior of superfluid 
""He, confined in aerogels, has not been observed. 
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FIG. 1: Relative critical temperature Tc(c)/Tc(0) (circles) compared to the mean-field approximation (dashed line). The 
dotted line corresponds to the random dilution case. 
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FIG. 2: Finite size scaling of maximum logarithmic derivatives (a) and helicity modulus (b), for the pure 3DXY and c = 
0.05, 0.1, 0.2 as indicated. Plots have been displaced vertically for clarity (vertical axes are in axbitray units). Corresponding 
FSS slopes are indicated on each graph. 
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volume fraction c, taken from FSS fits in Fig. 2 (black squares), and from FSS of 
maximum specific heat (black circles), using Josephson's relation u* = 2/{d+a/u). Results from Refs.P, Ql^li depicted here for 



FIG. 3: Exponents u and v 
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comparison (corresponding open symbols). Dotted lines show average values for c 7^0 as indicated 
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FIG. 4: Plot of n(r) versus r for structures within DLCA, as indicated, at c = 0.05 and 0.2. Arrows indicate crossover to 
homogeneous regimes in DLCA. Lattice sizes are L = 64, plots result from averaging over 32 realizations. The dashed line 
corresponds to the slope —l/v. 



